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Abstract. We report a new analytical method for solution of a wide class 
of second-order differential equations with eigenvalues replaced by arbitrary 
functions. Such classes of problems occur frequently in Quantum Mechanics 
and Optics. This approach is based on the extension of the previously reported 
differential transfer matrix method with modified basis functions. Applications 
of the method to boundary value and initial value problems, as well as several 
examples are illustrated. 



1. Introduction 

Analytical solution of differential equations provides insight to the behavior of 
solutions whenever they exist. Unfortunately, most practical physical problems 
are described by governing models which are normally solved by numerical tech- 
niques. Moreover, explicit solutions even to the simplest differential equations are 
rare. A wide class of physical problems are described by second-order differential 
equations, where the only known analytical method with explicit solution for this 
purpose is the approximate Wentzel-Kramers-Brillouin (WKB) method [HE]. An 
approximate method for solution of linear homogeneous differential equations with 
variable coefficients has been reported in [3] , which is based on a transformation into 
a Volterra or Fredholm integral equation. Also, a matrix method has been reported 
in [3] which is used to transform the Schrodinger equation into a ShabatZakharov 
system of second order, and then the solution is obtained by perturbation tech- 
nique. There are also a number of further existing analytical methods for solution 
of Ordinary Differential Equations which are categorized and reviewed by Polyanin 
and Zaitsev in [5|. 

Recently, we have introduced an analytical method which is capable of solving 
linear homogeneous Ordinary Differential Equations (ODEs) with variable coeffi- 
cients [6]. The method is based on the definition of jump transfer matrices and 
taking the differential limit. The approach reduces the nth-order differential equa- 
tion to a first order system of n linear differential equations. The full analytical 
solution is then found by the perturbation technique, which may be elegantly ex- 
pressed in terms of matrix exponentiation of the integral of a Kernel matrix. The 
important feature of this method is that it deals with the evolution of indepen- 
dent solutions, rather than its derivatives. The exact mathematical nature of this 
method has also been rigorously established [B]. 

This method emerged as an extension of Differential Transfer Matrix Method 
(DTMM), which was originally proposed in the context of optics and quantum 
mechanics [3 El [9] . DTMM is based on the modification of the standard transfer 
matrix method in optics [10] and quantum mechanics [11] , and employs exponential 
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basis functions. Through numerous examples, this method was shown to be simple, 
exact, and efficient, while reflecting the basic properties of the physical problem. 

However, the initial formulation of DTMM had difficulty when dealing with 
singularities. Such singularities arise in the domain of optics and quantum mechan- 
ics at the turning points of wavefunctions, where the behavior switches between 
oscillatory and decaying forms f2l. As a result, basic DTMM results suffer from nu- 
merical error when approaching physical singularities. To overcome this difficulty, 
the solution may be expressed by using Airy functions [5]. Further applications 
and extensions of DTMM have been reported by Mehrany et. al., which include 
employment of WKB basis |12i 113] for improvement of numerical accuracy, and 
numerical implementation of the method using Airy functions [M] for dealing with 
singularities. Also, the original DTMM has been used in [T5| with no essential 
change as reported in the main formulation [^. 

In this paper, we report a general formulation of DTMM for eigenvalue problems 
having second-order ODEs. The generalization is done by replacing the eigenvalue 
with an arbitrary function, which we here refer to as the eigenvalue function. We 
show that DTMM is still capable of delivering the analytical solution provided that 
proper basis functions are used. We also establish the mathematical validity of this 
method by presenting the fundamental theorem of DTMM. The major improvement 
of this work in contrast to the original formulation 6 is two-fold. Firstly, the 
solution to the extended eigenvalue problem evolves 'naturally' out of the relevant 
basis functions rather than exponentials. Secondly, envelope functions undergo 
minimal variations because most of the solution is encompassed in the extended 
bases. As a result, we may obtain a simple approximate solution for problems with 
slowly varying eigenvalue functions. 

2. Formulation 

Suppose that a linear homogeneous ordinary differential equation (ODE) is given 

as 



Mkf{x;k)=0 (1) 
where Hfe is a second-order linear operator with eigenvalue k which involves partial 
derivatives only with respect to x, and / (x; k) being its eigenfunction. It is the 
purpose of this formulation to find solutions of ([!]) when the eigenvalue k is replaced 
by an eigenvalue function such as k (x). Hence, the extended problem reads 

Lfe(,)/[x;fc(x)] =0 (2) 
Here, the new operator h^^) is obtained from Hfc by replacing the eigenvalue k 
with an eigenvalue function k{x), and d/dx with d/dx. In general, ^ admits a 
general solution having the form 

fix;k) = aA{x;k)+bB{x]k) (3) 

in which a and b are constants determined by initial or boundary conditions, and 
A (x] k) and B (x; k) are linearly independent solutions with non- vanishing Wron- 
skian determinant, i.e. 



W{x;k) = 



A{x;k) B{x;k) 
A^ (x; k) B^ {x; k) 



7^0 



(4) 
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Clearly, Ax (x; k) — dA {x; k)/dx and Bx [x; k) — dB {x; k)/dx. 

Now, we seek a solution to ([2]) having the extended form to ([3]) given by 

/ (x; k)=a {x) A [x\ k {x)] +b{x)B [x; k (x)] 
where a (x) and b (x) remain as unknown functions to be determined. 



(5) 



2.1. Transfer Matrix of Finite Jumps. In order to proceed with the formula- 
tion, we take a similar approach to the conventional Differential Transfer Matrix 
Method with exponential basis. For this reason, we first need to obtain the transfer 
matrix of finite jumps in the eigenvalue k. Suppose that the eigenvalue function is 
k (x) defined as 



k{x) 



fci, X > X 
k2, X < X 

with fci, k2, and X being constants. Solution to (H]) is readily given by ([3]) as 



(6) 



/ [x; k (x)] 



fi{x), x<X 
/2 (x) , x> X 
aiA (x; ki) + hiB (x; fci) , 
a2A (x; fc2) + b2B (x; fc2) , 
Analyticity of / [x; fc (x)] across requires that 



x<X 
x> X 



(7) 



/i (X) = f2 (X) 
fi{X)^f^{X) 
We hence arrive in the system of equations 



(8) 



a2 
b2 



(9) 



A{X;k2) B{X-k2) 
Ax{X-k2) Bx{X;k2) 

AiX;ki) B{X;k,) 1 f a, 
Ax{X;k^) Bx{X-M) \ I &i 
Since the Wronskian W is supposed to be non-zero by dH), then © can be solved 
to obtain 



02 

b2 



1 

W2 



A2 B2 ' 


-1 


" ^1 


Bi ' 


A'2 B'^ 






B[ 


B'2 ~B 


2 






Bi 


-A'2 M 








B[ 



ai 

bi 
ai 
bi 



Here, Ai and A'^ respectively denote A {X; fcj) and Ax {X: ki) 
Bi and B'i respectively denote B{X;ki) and Bx{X;ki), i 
W {X;k2)- Hence, we get 



(10) 



i = 1,2. Similarly, 
= 1,2. Also, W2 = 



{C2} = [Q'^'] {C,} = 



1^2 1^2 
yii 112 

I21 I22 



{Ci} 



(11) 



Here, [Q^^^] is referred to as the jump transfer matrix with the elements 



„1^2 

111 
1^2 

921 



W2 

A{A2-AiA'2 
W2 



q'2^' 



W2 

B[A2-BiA'2 
W2 



(12) 



and 



4 



SINA KHORASANI 



j,^ = l,2 (13) 

Properties of jump transfer matrices have been extensively discussed in earlier 
works. But for the sake of convenience we mention a few 



lni^2| _ Wi_ 
\^ \ ~ W2 

corresponding respectively to the determinant, combination, and inversion proper- 
ties. Moreover, we readily notice that = Wm/Wn- It is clear that a given 
transfer matrix [Q^^'^] is not invertible unless the Wronskian ^ does not vanish. 
The combination property explains how to obtain the total transfer matrix over a 
number of finite jumps, among which the eigenvalue function k (x) is constant. 



2.2. Differential Transfer Matrix. Now, we let k {x) be a smooth function of 
X. Within the infinitesimal neighborhood of any given point such as x = X, the 
eigenvalue function k (x) will undergo a first-order change from ki — k {X) to 
k2 = k{X + Ax). We may then define k2 = ki + Ak, where Ak represents a small 
change in the eigenvalue. If Ax is small, then we may neglect the variations of k (x) 
within [X, X + Ax]. The corresponding first order change in the vector {C} across 
X = X will be clearly given by 

{AC}«^([Qi-2]-[/]){Ci} (15) 

where the approximation becomes exact if we let Ax approach zero. Thereby, we 
get 

d{C{x)} :^[U{x)]{C{x)}dx (16) 
in which {C (x)} is the envelope vector function, and the Kernel matrix [U (x)] is 



We notice that in order to obtain the correct solution to the Kernel matrix 
[U (x)], one needs to make the replacements Gi — G, G[ — Gx, and G2 = G + 
GfeAfc, G2 — Gx + GxkAk, where represents either of A or B. Here, subscripts 
refer to partial derivatives in the sense that Gx = dG/dx, Gk — dG/dk, Gxk = 
d^G/ dxdk. After doing some algebra and simplification we get the complete form 
for the elements of the Kernel matrix [U (x)] as 

U2^^k'^^^-^^ U22^k'^-^^^^ ^^^^ 

Here, k' — dk {x)/dx and W = ABx — AxB. In the fully expanded form we have 



[U{x)] = 



Uii (x) Ui2 {x) 
U21 (x) U22 (x) 
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-1. (-) = #S I^^S^^ (-)] - »M»M} 

(.) = #i I M^M^ _ ^2^A fc (X)]} 



(19) 



where 



W{x) = A fc (.)] ^^MM _ dA[x,kix)] ^ ^ 
ox ox 
A general solution to ^ is given by ([S|) with 

{C'ix)} = [Uix)]{Cix)} (21) 

in which the elements of the Kernel matrix [U (x)] are given in ([T8|. Interestingly, 
(|2T]) allows an analytical solution through perturbation theory as [2 [6] 

{C {X2)} = {C (x^)} +J[U (yo)] {C iyo)}dyo 

XI 

X2 yi 

+ 1 IIU (yi)] [U (yo)] {C (yo)} dyodyi (22) 

3^2 2/2 yi 

+ / / / [C^ (y2)] [f^ (yi)] [U (yo)] {C (yo)} dyodyidy2 + ■■■ 

Xi Xi Xi 

Often ([22]) is written symbolically as 



{C{x2)}^Texp^J^[U{x)]dxj{Cixi)} 
= {C(xi)} + /[C/(yo)]{C(yo)}dyo 

X2 X2 

+i / / T [C/ (yi)] [C/ (yo)] {C (yo)} dyodyi 



(23) 



Xi Xi 
X2 X2 X2 



LJ J JT[U (y^)] [U (yi)] [C/ (yo)] {C (yo)} dyodyidys 



' 3! 

xi xi xi 

= [Q^i^"^]{C(a;i)} 

where [Q^"^"*^^] is the transfer matrix from xi to X2 and exp(-) being the matrix 
exponentiation 

oo ^ 

exp[M] = [/] + ^-[M]" (24) 

n—l 

Furthermore, T is the Dyson's ordering operator defined as 



T[a(«)l[aW|.{P(°)l[^(„')j; ,25) 

There are few known sufficient conditions [51 116] , which rarely happen and under 
which T may be dropped exactly to reach 
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{C(x2)} = exp 




{C(xi)} 



(26) 



The most important sufficient condition includes the case when the Kernel matrix 
commutes with itself as \U (a)] \U (6)] = [U (6)] \U (a)] for all given a and h. This 
is better known as the Lappo-Danilevskii [TT] criterion, and is also generalized by 
Fedorov [18]. Evidently, this condition applies to all constant as well as diagonal 
Kernel matrices. It is therefore a challenge to find construct the proper extended 
eigenvalue equation in such a way that the corresponding Kernel matrix meets any 
of these sufficient criteria. If possible, then exact and explicit closed form solutions 
to ([2]) arc found analytically by using ([26]) instead of (|23|). 

While in general, ([26]) is merely an approximation to the exact solution ([23]) . 
nevertheless, in any case the determinant and trace of [Q^^"*''^^] would be preserved 
exactly, meaning that the eigenvalues of the transfer matrix [Q^^"*^^] would not 
be affected at least. This allows us to formulate a very convenient approximate 
numerical solution to the extended problem ([2]). 



2.3. Properties of DifTerential Transfer Matrix. Properties of the transfer 
matrix [Q^^'^^'^] as defined in ([23l) are very much similar to those of the jump 
transfer matrix (|14|) . given by 



JQX1^X2] ^ [Q^3^X2] 



(27) 



Properties of the transfer matrix in which Wi = W [xi; h] — W [xi; k (xi)] , i — 
1,2. The first property, the identity property readily follows by definition in (1231) . 
The determinant property can be observed by noting that Dyson's operator has no 
effect on the determinant, and thus 



Texp<^ / [U{x)]dx 



exp \ J [U (x)] dx 



exp tr < J [J7 {x)] dx>\— exp J [uu {x) + U22 [x)] dx 



(28) 



Now from (fT8|) we get 



7 [^11 (x) + U22 {x)] dx = 7 \^lz^f: k'dx = - / i In W^dx 

fc2 (29) 
= - / ^ In Wdk In fci - In /c2 

fci 

Inserting this result in (|28p gives the determinant property in (|28p . Combination 
and inversion properties also follow the definition 
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2.4. Fundamental Theorem. In this section, we rigorously show that the differ- 
ential transfer matrix method leads to an exact solution of the differential equation 
To show this, we start by proving the so-called Derivative Lemma. 

Lemma 2.4.1. The total derivative of function is given by 



9" 9" 
— / [x; k (x)] = a {x) —A [x; k {x)] + b (x) —B [x; k {x)] , < n < 2 (30) 

Proof. Proof follows by direct substitution. The case of n = is trivial, and we 
first prove the validity of (|30|) for n = 1. Direct differentiation of ([5]) by chain rule 
gives 

= a' A + b'B + a [A., + fc'Afe) +b{B^ + k'B^) (31) 

ax 

where / = aA + bB. But we already have obtained the derivatives a' and b' from 
(HU and (UHl) as 



a' = Una + ui2b = ^ [{A^kB - A^B^) a + {B^kB - BkB^) b] 
b' = U2ia + U22b = ^ [{AkA^ ~ A^kA) a + {BkA^ ~ B^^A) b] 

in which the elements Uij,i,j — 1, 2 of the Kernel matrix are given in After 
some minor algebra we get 

^/ = aA, + bB, (33) 

To show the correctness of ((30|) for n = 2, we take the derivative again with respect 
to X from both sides of ([30]) . thus giving 



(34) 



£,f^£(aA,+bB,) 

= a' A, + b'B, + a±A, + b£B, 

= [una + ui2&) A, + (u2ia + U22b) B, 

+a {A,, + k'A,k) + b [B,, + k'B,k) 



Further substitution of ([5^ in ([M]) gives 



■&f = abJ'-ba^ [i^^^kB - AkB,) a + {B,kB - B^B,) b] A, 

+ abXba^ MkA, - A,kA) a + {BkA, - B,kA) b] B, (35) 
+a {A,, + k'A,k) + b [B,, + k'B,k) 

Here, we have used the definition of the Wronskian W — AB, — BA, from Q. 

Eventually, after some algebra we arrive at the final result 



d2 

This completes the proof 



-^f = a A,, + bB,, (36) 



□ 



We are now in a position to present the Fundamental Theorem of Differential 
Transfer Matrix Method as follows. 

Theorem 2.4.2. The solution to 0) having the form (0) with Ii21\} is exact. 
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Proof. To show the vahdity of the statement, we start by plugging in ^ directly 
into ©. 



(37) 



Lfc/ [x; k] ^ Lfe {aA [x; k] + bB [x; k]} 
= oRkA [x; k] + bUkB [x; k] 

where Hfc and Lfc are respectively defined in (IT|) and ©. Here, the explicit de- 
pendence of the operator Lfc and functions a, b, and A: on a; is not shown. But by 
assumption we have M^A [x; k] — MkB [x; k] — and hence the assertion. □ 

All remains is to force the initial or boundary conditions, and we here mention 
the treatment of both types of conditions. 

2.5. Initial Conditions. Without loss of generality, we may assume that initial 
conditions are known at some point like x = c. Suppose that / (c) and /' (c) are 
known. Then by derivative lemma ((30|) we have 



/' (x) = A^ [x; k (x)] a {x) + B^ [x; k {x)] b {x) 



Therefore we get 



/(c) 
fie) 



A [c; k (c)] B [c; k (c)] 
A,[c;k{c)] B,[c;kic)] 



Thus we may obtain the initial vector {C (c)} as 



{C{c)}^ 
with the solution 

{C{c)} 



A [c; k (c)] B [c; k (c)] 
A,[c;k{c)] B,[c;k{c)] 



a{c) 
6(c) 



/(c) 
/'(c) 



1 



B, [c; k (c)] -B [c; k (c)] 
-A^ [c; k (c)] A [c; k (c)] 



W [c; k (c)] 

Now by ((23)) we have the full solution to the envelope vector as 



/(c) 
/'(c) 



(38) 
(39) 

(40) 
(41) 



(42) 



{C(a;)} = Tcxp|/[{/(y)]dy|{C7(c)} 

«exp|/[C/(y)]dz/|{C(c)} 
This allows us to obtain the final solution to the initial value problem ([2]) as 

f{x)=[A [x; k (x)] B [x; k (x)] ] {C (x)} (43) 

2.6. Boundary Conditions. Without loss of generality, we may assume that 
boundary conditions are known at two points like x = ci,C2. Suppose that / (ci) 
and / (02) are known. Then by ^ we have 



/(ci) = 

/(C2) = 

But from (l23l) we have 



A[ci;k{ci)] B[ci;k(ci)] 
A[c2;kic2)] B[c2;k{c2)] 



{C{ci)} 

{C{C2)} 



(44) 



{C(c2)} = [Q^^-^^]{C(ci)} (45) 
Thus by combining and pS)) we may obtain the system of equations 
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/(ci) 

/(C2) , 

A[ci;k{ci)] B[ci;k{ci)] 
qiiA [c2; k (C2)] + q2iB [02; fc (02)] QuA [02; k (02)] + 922-8 [02; k (02)] 

(46) 

with the solution 



{C{c^)} 



A[ci;fc(ci)] B[ci;fc(ci)] ' ' 



X 



(jiiA [c2; k (C2)] + (J21-B [c2; k (02)] 512^ [c2; fc (02)] + (722B [02; (02) 
I /(ci) 

/(C2) 

{C(c2)} = [Q^^-^^]{C(ci)} 

(47) 

The rest is the same and similar to (H^ 

{C (x)} = Texp I / [f/ (y)] dyj {C (ci)} = Texp | - / [[/ (y)] dyj {C (02)} 
« exp |/ [[/ (j/)] dyj {C (ci)} « exp |- / [(7 (y)] dj/j {C (02)} 

(48) 

3. Examples 

In this section, we present some application examples describing the details of 
the method. 

3.1. Wave Equation. Numerous physical problems are described via the simple 
second-order equation 

V-" (a;) + k^ (x) ip{x)^0 (49) 

The equation (|49p is known to have no explicit solution for arbitrary eigenvalue 
function k (x) , which in the literature is usually referred to as the wavenumber 
function. The only existing analytical solution to (|49l) is the very well-known WKB 
approximation. Here, the corresponding operators read 

In case of constant wavenumber, any solution of (|49|) is given by linear combinations 
of exponential functions exp(±ifca;), and therefore the eigenfunctions are readily 
found to be 

A [x; k (x)] = exp [—ixk (x)] B [x; k (x)] = exp [+ixk (x)] (51) 

with the Wronskian W \ 2ik {x) . To find a solution to (|49|) it is therefore 

sufficient to find the Kernel matrix \U {x)] , which is 
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[U{x)] 



(52) 



fc'(x) r -l + 2ixk(x) exp[+2ixk{x)] 
2k (x) exp [—2ixk (x)] —l — 2ixk{x) 
Numerical and analytical solutions of (|49p using the Kernel matrix ([5^ has been 
extensively studied in the literature through past years. 



3.2. Airy's Equation. One of the known problems with the above differential 
transfer matrix solution of p9|) is occurrence of singularities at which fc^ (x) changes 
sign and the Wronskian vanishes, better known as returning points. The singularity 
is clear in ([5^ as the denominator. Even the approximate WKB solution fails near 
returning points. It is a common practice however to expand the solution near the 
returning points via Airy functions. Without loss of generality we may assume that 
the singularity is located at a; = 0. In this case, we present a slightly modified form 
of (mi) as 



i;" (x) - k^ (x) xij) {x) = 
with the Airy wavenumber satisfying k (0) ^ 0. Hence, 



(53) 



- k^ {x)x 



(54) 



0)2 .. 

In case of constant wavenumber, any solution of ()53|) is given by linear combi- 
nations of Airy functions Ai (kx) and Bi (kx), and thus the eigenfunctions take the 

form 



A [x; k (x)] = Ai [xk (x)] B [x; k [x)] = Bi [xk (x)] (55) 

having the constant Wronskian W [x; k (x)] = I/tt. This choice of basis, after some 
simplification, leads to the following Kernel matrix elements 

uii = +TTk' {fc^x^Ai (fcx) Bi (fcx) + A*' {kx) [Bi (kx) - kxBi' (kx)]} 
ui2 = irk' [Qi' {kx) [Bi (fcx) - xfcB«' (fcx)] + fc^x^Bi^ (fcx)} , . 

= Trfc' {Ai' (kx) [xkAi' (kx) - Ai (fcx)] - fc^x^Ai^ (fcx)} ^ ' 

= -Trfc' {fc^x^Ai (fcx) Bi (fcx) + Bi' (fcx) [Ai (fcx) - fcxAi' (fcx)]} 

which is clearly non-singular at x = 0. Here, dependence of the eigenvalue function 
fc (x) on X is not displayed for the sake of convenience. 

3.3. Bessel's Equation. In the domain of fiber optics having cylindrical symmetry 
and after proper transformations, the radial component of the wave equation takes 
the form 

xV (x) + xij' (x) + [fc2 (x) x^ - 1^2] ij{x) = (57) 

Hence 



Hfc=x2|^+x^ + (fcV-^.2) L,(,)=x2-^+x-^-f[fc2(x)x2-t.2i (58) 
ax^ ox rfx^ ax 

In case of constant wavenumber, any solution of (j57p is given by any linear 
combinations of Bessel and Neumann functions (fcx) and N^(kx), or Hankel 
functions Hi^'' (kx) and Hi'^^ (kx). For the first pair the eigenfunctions take the 
form 



GENERALIZED DIFFERENTIAL TRANSFER MATRICES 



11 



A [x] k (x)] = J„ [xk (x)] B [x; k (x)] = [xk (x)] (59) 

having the Wronskian 14^ [x; fc (a;)] ~ 2/t:x. The elements of the corresponding 
Kernel matrix are 

uii - +^ {[kx [J,_2 {kx) - 2 J, {kx) + J,+2 {kx)] + 2 [J,_i (kx) - J,+i (fca;)]] N, (kx) 

-kx [Ju-i {kx) - Ju+i {kx)] [N^^i {kx) - N^+i {kx)]} 
ui2 = {[kx [iV,_2 {kx) - 2N, {kx) + N,+2 {kx)] + 2 [iV,_i {kx) - iV,+i (fcx)]] iV, (fcx) 

+fcx [iV^_i (fcx) - 7V^+i (fcx)]^| 
U21 = +^ {[fca; [J^_2 {kx) - 2J^ {kx) + {kx)] + 2 [J^_i {kx) - J^+i (fcx)]] {kx) 

+kx [Ju-i {kx) - Jy+i (fca;)]^| 

U22 = {[kx [iV,_2 {kx) - 2N, {kx) + N,+2 {kx)] + 2 [iV,_i {kx) ~ iV,+i {kx)]] J, {kx) 

-kx [J^^i {kx) - J^+i {kx)] [N^-i {kx) - N^+i {kx)]} 

(60) 

As it can be seen, all elements of the Kernel matrix vanish at x = 0; this ensures 
cylindrical symmetry of the radial function. 

It should be also noticed that since Bessel functions with non-integer order are 
well-defined [19] 

^-(•-)-E ^,r(L + . + i) (2) (^1) 

One could have also replaced the eigenvalue function with the order as 

x'^ip" {x) + xil)' {x) + [a;2 - k^ {x)] tj} {x) = (62) 

Hence 



Again since Jq, (x) and J_ct (a^) for non-integer a are linearly independent, the 
natural choice of eigenfunctions based on Bessel functions would be either 



A [x] k {x)] = J+k(x) {x) B [x] k {x)] = J^k(x) {x) (64) 



or 



A [x] k {x)] = Jk[x) {x) B [x; k {x)] = N^^) {x) (65) 

The latter pair has the same Wronskian of W[x]k{x)] = 2/-kx. Hence, the 
elements of the Kernel matrix are 



.11 (X) = ^ {^^^.(.) (X) - 'J^'J^^] 

"12 yx) — 2 I dxdk ^^k(x) yx) gj, j 

_ Tixk'jx) f 8Jfc(^)(a) 3^Jfc(^)(a:) ^ , .\ 

"2U-i; — 2 \ dk dx dxdk 'Jk{x)(X)j 

U22 {X) - — 2 — \ — 5^ dk •-'H^) {^) dxdk — / 



(66) 
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3.4. Euler-Cauchy Equation. The Euler-Cauchy equation reads 



X^V" (x) + Xl/y (x) - {x) %P {x) 

Therefore, the corresponding operators are given by 



Wk^x"- 



x^-k^ 
ox 



Li-i-T-i = X 



X— k^ (x) 



(67) 



(68) 



dx^ ' dx dx'^ ' dx 

When k (x) is a constant, solutions of (p7|) are given by linear combinations of x^^ 
and a;"'^. Thus, the extended basis functions are 



A [x; k {x)] ^ x+'^'f^) B [x; k {x)] = x^'^^^^ (69) 

with the Wronskian W[x;k{x)] = ~2k (x) /x, which vanishes at k{x) — 0. The 
Kernel matrix simphfies to the convenient form 



[Uix)] 



k' (x) 
'2k{x) 



l + 2fc(a;)lnx 

_^2fc(x) l-2fc(a;)lna; 



(70) 



3.5. Approximate Solution. The exact matrix exponential in can be evalu- 
ated if off-diagonal elements of the Kernel matrix could be dropped. This has been 
previously shown to lead to the well-known WKB solution [H [TJ |H] . Under such 
conditions, the approximate differential transfer matrix takes the form 



[gxi^x2] ^Texp<^ / [U{x)]dx I wTexp<^ / 



Mil {x) 

M22 (a;) 



dx 



exp 



X2 

J Mil (x) dx 

xi 









exp 



/ M22 (x) dx 



(71) 



3.6. Periodic Perturbations. There is a great deal of simplification possible, 
when the eigenvalue function k (x) = k (x + L) is periodic for some L > 0. Let 
Tl be the translation operator, defined as TLh{x) = h{x + L), for all arbitrary 
functions h{x). Hence, we readily have the commutation property [L/;(j,),Tl] = 
0, and hence these two operators share identical eigenfunctions. Since L^^^,) is 
supposed to be linear, Bloch-Floquet theorem |9j applies and then any solution will 
take the form of Bloch eigenfunctions 



Ti/ [x; k {x)] ^ f[x + L;k{x + L)]^ exp {~jkL) f [x; k (x)] (72) 

in which the complex number k is being referred to as the Bloch number. This 
shows that / [x; k (x)] is an eigenfunction of the translation operator with the 
eigenvalue exp (—juL). Based on (I72p . we furthermore have 

/ (x) = exp {-jnx) (x) (73) 
with the envelope function satisfying 



TiffK {x) ^ 5k {x) 



(74) 
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Under such circumstances, it is easy to obtain the characteristic equation of 
eigenfunctions. From ([5]) and ([72l) we have 



TlJ [x; k {x)] ^a{x + L)A[x + L;k (x)] + b {x + L) B [x + L; k {x)] 
= exp {—jnL) {a (x) A [x; k {x)] +b{x) B [x; k (x)]} 



(75) 



By taking differentiating with respect to x from l|72p and derivative lemma pop we 
also get 



(76) 



Ti/' [x; k (x)] = exp (-jkL) f [x; k [x)] 

= exp {—JkL) {a (x) Ax [x; k {x)] + h (x) B^ [x; k (x)]} 

We furthermore notice that ((75|) and ([75]) are actually identities which hold for 
all X. These two can be combined to get the system of equations 



exp {-jnL) 
which allows the solution 



A[x + L;k{x)] B[x + L;k{x)] 
Ax [x + L;k (x)] Bx [x + L;k (x) 

A[x;k{x)] B[x;k{x)] 
Ax [x; k (x)] Bx [x; k (x)] 



{C{x + L)} = 
{C{x)] 



(77) 



ilL yx)l — ^[x+L-Mx)] 



Bx[x + L;k{x)] ~B[x + L;k{x)] 
-Ax [x + L\k [x)] A[x + L]k {x)] 
x; k (x)] B [x; k (x)] 
Ax [x; k (x)] Bx [x; k (x)] 
e^p{^jKL)[V] {C{x)} 



A 



{C{x)} 



(78) 



But from ^ 



Tl{C{x)}^[Q^^^+^] {C{x)} 
Simultaneous satisfaction of ([77)1 and (1751) requires that 



(79) 



\exp{-jKL)[V]~[Q]\^Q (80) 
In other words, we should have 

«=^ln{eig[P]} (81) 

in which [P] = [V]^^ [Q]. Here, constancy of the Bloch wavenumber k is guaranteed 
by Bloch-Floquet theorem, and is rigorously shown to hold for the example of 
extended wave equation (|49p with periodic wavenumber in [9]. 



4. Conclusions 

In this paper, we presented a new analytical solution obtained by differential 
transfer matrix method to a wide class of second-order linear differential equa- 
tions, which are extended from eigenvalue problems by replacing the eigenvalue 
with an arbitrary eigenvalue function. We presented the details of the method and 
a fundamental theorem to rigorously establish the mathematical formulation. Few 
examples were also described. 
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